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o , 

We calculate the electron-electron interaction induced energy-dependent inelastic carrier relax- 
, ation rate in doped semiconductor coupled double quantum well nanostructures within the two 

subband approximation at zero temperature. In particular, we calculate, using many-body theory, 
the imaginary part of the full self-energy matrix by expanding in the dynamically RPA screened 
Coulomb interaction, obtaining the intrasubband and intersubband electron relaxation rates in the 
ground and excited subbands as a function of electron energy. We separate out the single particle 
and the collective excitation contributions, and comment on the effects of structural asymmetry in 
the quantum well on the relaxation rate. Effects of dynamical screening and Fermi statistics are 
automatically included in our many body formalism rather than being incorporated in an ad-hoc 
, manner as one must do in the Boltzman theory. 
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I. INTRODUCTION 



Electron-electron interaction induced carrier relaxation is an important inelastic scattering process in low- 
dimensional semiconductors nanostructures. It is often (particularly in situations where LO phonon emission is 
energetically prohibited because the excited electrons do not have enough energy) the most dominant relaxation pro- 
cess in semiconductor quantum wells and wires, and is therefore of considerable fundamental and practical importance. 
Band gap engineering has led to the possibility of fabricating tunable far infrared quantum well cascade lasers (QCL) 
and efficient quantum well infrared photodetectors (QWIP) where inelastic carrier relaxation via electron-electron in- 
teraction is a crucial (perhaps even decisive) process in determining device operation and feasibility [jjj . For QCL and 
QWIP operations it is the intersubband inelastic relaxation which turns out to be the primary rate limiting scattering 
process. For other proposed devices, such as the planar hot electron transistors or related 2D high speed devices, 
intrasubband relaxation is the important process. A thorough quantitative understanding of intra- and inter-subband 
relaxation due to electron-electron interaction is therefore important for the successful realization of these devices. 
In addition to this practical technological motivation arising from QCL, QWIP, and other proposed band-engineered 
quantum well devices, there is also obvious fundamental reason for studying inelastic Coulomb scattering in 2D quan- 
tum well systems. Inelastic electron-electron scattering determines the 2D quasiparticle spectral width, as determined 
for example in tunneling measurements, through the imaginary part of the electron self-energy function 
■ In this article we use a many-body approach in calculating the inelastic relaxation rate of 2D electrons confined in 
GaAs-Al^Gai-a; As semiconductor quantum well structures. Our work is a multisubband generalization of the earlier 
work H by Jalabert and Das Sarma who considered only intrasubband relaxation within a single subband model. We 
' consider both intra and intersubband relaxation in the two lowest subbands, and consider both single well and coupled 
^ . double well structures. An additional important issue addressed in our work is the effect of structural asymmetry 
in the quantum well on the relaxation rate. This is in fact a potentially significant factor in the fabrication of QCL 
and QWIP structures since asymmetry could lead to the opening of new electron-electron interaction channels in the 
inelastic intersubband relaxation as we discuss below in this article. 

The central quantity we calculate in this work, within the leading order dynamically screened Coulomb interaction 
expansion (the so-called GW approximation in the multisubband situation), is the imaginary part of electronic on-shell 
self-energy matrix, M, in the quantum well subband index etc.). The subband self-energy in the multisubband 
situation is, in general, off-diagonal reflecting the breaking of the translational invariance along the growth (2) direction 
(we take the x-y plane to be the 2D plane with all wave vectors in this paper being 2D wave vectors in the x-y plane). 
The off-diagonal self-energy, Im(-My), incorporates in an intrinsic many-body manner the possibility of electron- 
electron interaction induced intersubband scattering (both virtual and real) of carriers. We believe that in the 
doped situation of our interest, where the quantum well subbands are occupied by many electrons, the many-body 
self-energy approach is the only reasonable technique in calculating the inelastic carrier relaxation rate in contrast 
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to the Boltzman equation approach where the scattering rates usually calculated using Fermi's golden rule. This is 
because the dynamical screening inherent in the many electron system, which affects the calculated inelastic scattering 
rates in profound and highly nontrivial way, is automatically incorporated in our many-body GW expansion whereas 
inclusion of dynamical screening in Fermi's golden rule type formula is problematic and can only be done heuristically 
by replacing the bare interaction by a screened interaction in an ad-hoc manner. 

Our theory, as mentioned above, is based on the so-called GW self-energy approximation where the electron self- 
energy M is obtained in a leading order expansion of the dynamically screened Coulomb interaction W = V s , where 
the superscript s denotes dynamical screening of the bare electron-electron interaction matrix V in the multisubband 
situation. We use the RPA to obtain the dynamical screened interaction V s , i.e. V s = e~ l V , with e = 1 — VII, where 
II is the leading-order (i.e. noninteracting) electron polarizability matrix. We also approximate the electron Green's 
function G by the noninteracting Green's function G° making our formal expression for the self-energy matrix to be 



M ~ J G U V S , (1) 

where the integral involves integrating over all internal momentum and energy variables as well as summing over all 
internal subband indices (and spin). Putting the subband (matrix) indices explicitly in Eq. (|l|), we get 

ImMtf=Im£; j G° lm V? lmj . (2) 

Im 

We note, however, that G°, being the noninteracting Green's function, is necessarily diagonal in subband indices (i.e. 
an electron cannot suffer an intersubband scattering in the absence of interaction); 

Ghn ~ G u 8i m . (3) 

Then, Eq. (2) becomes 

haM ij= J2 [lm[G° u Vi Uj \, (4) 
i J 

with 

Equations (|J) and (^) are the central formal equations we use in our theory to obtain the inelastic relaxation time r, 
remembering that the scattering rate T and the relaxation time r are connected by 

"5P (6) 



where, 



r=|ImM|. (7) 



In Eqs. (0) and (@), F = |ImM| is calculated on-shell, i.e., the quasiparticle self-energy defines T. To demonstrate 
how Eq. (|J) may, in principle, differ from the simplistic Fermi's golden rule approach we consider the specific two 
subband model of interest to us in this paper. Then i, j, I = 1, 2 with only subband 1, the ground subband, and the 
subband 2, the first excited subband being considered in the theory assuming all other subbands to be substantially 
higher in energy making negligible contributions to the self-energies of the lowest two subbands. We also assume the 
square well structure to be symmetric, so that parity is a good quantum number in the problem which makes all 
" off-diagonal" interaction matrix elements vanish |^| by virtue of parity conservation with the only non-zero elements 
of V s being Vf nl , V 2 2 2 2i V 1 s 212 , and V{ 122 (note that V 1221 = V 2 S 112 and Vf 122 — V 2211 by symmetry). In this situation 
Eq. ([|) implies that 

lmM 12 = ImM 2 i = 0, (8) 

and, 

ImAf 22 ~ J Im [G° n V 2 s n2 + G° 22 V 2 S 222 ] . (9) 
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We note that the dynamically screened interaction matrix element V{ 122 is not explicitly present in Eq. (|9j). On the 
other hand, a naive Fermi's golden rule approach will explicitly include such a V{ 122 term, as indeed has been done in 
the literature ||, because it seems to arise from the direct Coulomb interaction V1122 between an electron in subband 
1 and an electron in subband 2 without any intersubband scattering. We mention, however, that dynamical screening 
of V2222 produces an effective V22 11 term in our theory since dynamical screening proceeds through virtual creation 
of electron-hole pairs. 

We have calculated the energy dependent inelastic relaxation rate at T — for a two-subband (1 and 2) GaAs- 
Al^Gai-zAs quantum well system with a total electron density N e — 2 x 10 11 cm~ 2 for the following five distinct 
situations, (i) Two coupled symmetric quantum wells of width 150 A each with interwell tunneling induced by a 
tunnelling barrier of height 228 meV and width 30 A. Here the lowest two subbands are the so-called symmetric 
(bonding) and antisymmetric (antibonding) levels with energies E-y — 15.35 meV and E 2 — 17.03 meV, respectively. 
The third level E 3 = 60.53 meV is sufficiently high to be ignored {E Fl = E F - E x = 4.28 meV; E F2 = E F - E 2 = 
2.61 meV), with both subbands 1 and 2 occupied by carriers. These results are presented in Section III. A below, (ii) 
Two coupled asymmetric quantum wells with interwell tunneling (the same as in (i) above), with one well of width 
150 A and the other of width 140 A leading to E x = 15.93 meV and E 2 = 18.55 meV (E F1 = E F — E\ = 4.75 
meV; E F2 = E F — E 2 = 2.13 meV). Again, the next excited subband E% — 62.86 meV is high enough to be ignored. 
These results are presented in Section III. A below, (iii) Two coupled identical symmetric quantum wells of width 
150 A each with no interwell tunneling (i.e. the interwell barrier is taken to be infinity) and with a barrier width 
of 30 A. Here, Ei = E 2 = 23.87 meV (this degeneracy arises because the two wells are identical and there is no 
tunneling), E F i = E F2 = E F — E\ = E F — E 2 = 3.44 meV, and the next subband E3 — 96 meV is sufficiently high in 
energy to be neglected. These results are presented in Section III.B below, (iv) The same as in the last case with no 
interwell tunneling but an asymmetric situation with the two wells being different. One with a width of 150 A and 
the other a width of 142.4 A so that the subband Fermi energies Epi = E F — E\ = 4.75 meV and E F2 — E F — E 2 — 
2.13 meV, which are the same as in (ii) above. In this situation E\ = 23.87 meV, E 2 = 26.49 meV (again E3 can 
be neglected). Furthermore, to keep the average distance between the two electron layers the same as in (ii), we 
choose a barrier of width 28.8 A. These results are presented in Section III.B below in comparison with those in (ii). 
(v) A single symmetric quantum well of width 300 A and a barrier height 228 meV, which leads to the lowest two 
subbands at E\ = 4.88 meV, E 2 — 19.51 meV, and the Fermi energy E F i = E F — E\ = 6.88 meV (with E F < E 2 , 
so that the second subband is empty). In this situation, the next excited subband, £3 = 43.74 meV, is high enough 
in energy to be neglected. These results are present in Section III.C below. Our reason for studying the five different 
classes of systems described above is that we are interested in understanding the effects of interwell tunneling and 
structural asymmetry on the electron relaxation rate. In particular, asymmetry breaks parity conservation making 
the off-diagonal matrix elements of Coulomb interaction (e.g. V1112, Vm, Vmij V2111, V2221, V2212, ^2122, V222 all 
of which are zero in the symmetric situation) nonzero, leading to new inelastic relaxation channels not present in 
symmetric structures. For the sake of brevity we present results for a single representative carrier density and well 
parameters in each of the five cases. Our theory could be easily generalized to obtain finite temperature relaxation 
rates. Note that our goal here is to provide a qualitative understanding of how various physical parameters affect 
Coulomb scattering rates in 2D quantum wells. 

The plan of this article is the following. In Sec. II we present a brief theory with working formulae; in Sec. Ill we 
provide our numerical results and discussions; we conclude in Sec. IV with a summary. 

II. THEORY 

Our basic theory is outlined in the Introduction, where the formal expression for the self-energies to be calculated 
were given. Our central GW-RPA expression |?J for the self-energy can be explicitly written out by using the 
noninteracting subband Green's function 

G%{lo, k) = % [u - Ei(k) + Ep}- 1 , (10) 

where ui is a complex frequency (h = 1), and Ei(k) = Ei + k 2 /2m* is the noninteracting subband one-electron energy 
dispersion. Using Eq. ( |To| ) in Eqs. (^) and (|^), and carrying out the internal frequency integration, and taking the 
imaginary part after the on-shell analytic continuation, we get 

Im M l3 (k) = JL J2jd\lm [Vfa (q, £* (k + q) - & (k)] x 

{9 (&(k) - &(k + q)) - 9 (-£,(k + q))} , (11) 
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where the on-shell subband energy £j(k) is given by 

= Ek(k) - E F , (12) 

and 6 (x) = (1), for a; < (>) 0, is the usual Heaviside theta function. The dynamically screened Coulomb interaction 
is given by (see Eq. ||) 

V£ lm = (e-'V)^, (13) 
with the multisubband RPA approximation j?| defined by the dielectric matrix 

Sijlm = (1 - % m n°J, (14) 

where Viji m is the bare Coulomb interaction matrix element in the subband representation, and II?., the noninteracting 
polarizability, is given by 

U ^)-- 2 J (2^0,-^ + 0) + ^)' (15) 

where /,-(k) is the Fermi distribution function in the ith subband. In this paper, we take the impurity scattering 
induced background broadening 7 as being a small phenomenological damping parameter which equivalent to be 
working in the clean limit. We are therefore restricting ourselves to high mobility quantum wells with small impurity 
scattering induced level broadening. 

Using Eqs. (p~l|-p~5|) it is straightforward to calculate the imaginary part of the on-shell self-energy. For the sake 
of completeness, we show below the detailed expressions for Im Mij in the GW approximation for the two-subband 
model: 

ImMn(fc) = onn(fc) + 01221 (ft), (16) 
ImMi 2 (fc) = 01112 (fc) + 0-1222 (ft), (17) 
ImM 2 i (ft) = 02111 (ft) + 02221 (ft), (18) 



and 



Here. 



ImM 22 (ft) = 02ii 2 (fc) + 02222 (ft). (19) 

CTllll(fc)= (2^ / d2 q{ Im ^im(q^)]^(-^)- e (-ei(k + q))]}, (20) 
Oi22l(fc) = T^2 J d 2 q{Im[y i s 221 (q,A + Wo )][0(-A-c o )-0(-6(k + q))]}, (21) 

oni2(ft) = ^2 J d 2 q{Im[V7 U2 (q,A)] [9 {-A) - 9 (-£(k + q))]} , 
01222^) = ^ J d 2 q{Im[V 1 a aa2 (q,A + wo)] [9 {-A - uj ) - 6 (-&(k + q))]} , 
o 2 lll(ft) = 7^T2 J d 2 q{Im[y 2 s nl (q,A-^o)] [9 {-A + UJ ) - 9 (-&(k + q))]} , 



(22) 



(23) 
(24) 



o 2 22l(ft) = 7^2 / d 2 q{Im[U 2 s 221 (q,A)][0(-A)-0(-e 2 (k + q))]}, (25) 
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oaii2(fe) = ^2 / rf 2 q{im[^ii 2 (q^-^o)][e(-A + ^ o )-0(-a(k + q))]}, (26) 

and 

(72222(^ = ^2 / rf 2 q{Im[K, s 222 (q,A)][0(- J 4)-0(-e 2 (k + q))]}, (27) 

where cl>o = i? 2 — _Ei is the subband energy difference and A = A(q, k) = (2kq cos r\ + q 2 )/2m* with r\ being the angle 
between k and q; and m* = 0.07m e being the GaAs conduction band electron effective mass. Now, we define the 
total inelastic Coulomb scattering rate for an electron with wavevector k ( i.e. an energy of k 2 /2m* with respect to 
the subband bottom) in the subband 1 and 2 as 

a x (k) = ImMn(fc) +ImAf 12 (fc) (28) 

and 

o-a(fc) = fmM 21 (fc) + Im M 22 {k). (29) 

It is important to realize that the screened potentials V^ lm for j ^ I do not appear in Eqs. ( p0| - p7j ) and, consequently, 
do not explicitly contribute to scattering rate. They are implicitly induced in the theory thr ough d ynamical screening 
H as discussed before. Furthermore, all screened interactions V^ lm involved in the Eqs. ( p0| - |27j ) are obtained from 
the relation between the bare electron-electron potential 

V ijlm (q) = — fdz [ dz'Mz)Mz)e- 9 \ Z - Z '\<t>l(z')<t>rn(z') 
?£o J J 

and the inverse matrix of the dynamical dielectric function eiji rn {q,uj) (see Eqs. ( p"3| ) and (|l4|), where the in- 
dices i,j,Z,and m = 1,2). These bare Coulomb potentials Viji m (q) are calculated here by using both the one- 
electron wavefunction <pi(z) and the subband energy Ei obtained through the numerical solution of the Schrodinger 
equation in the z direction for the specific quantum well confinement potential. Furthermore, the potential 
Vijim(l) can be separated into intra- and intersubband terms, and understood as follows: (i) Intralayer (intra- 
subband) interactions Vim(g) = Va, V 2222 (g) = Vb, and Vii 22 (g) = V 22 n(g) = Vc representing those scatter- 
ing events which the electrons remain in their original well (subband); (ii) Interlayer (intersubband) interactions 
^1212(9) = ^2121(9) = Vi22i(<z) = ^2112(9) = Vd representing scattering in which both electrons change their well 
(subband) indices; and (iii) Intra-interwell (subband) interactions Viii2(<?) = ^1121(9) = Vi2ii(g) = V2in(g) = Vj and 
^2212(9) = V222i(<z) = ^1222 (q) = ^2122(9) = Vh indicating the scattering in which only one of the electrons suffers 
the interwell (intersubband) transition. 

For each wave- vector k, the two dimensional q integrals in Eqs. d20|-p7|) are performed within the planes determined 
through the variables q and A in the screened interactions V^; m . The integration domains of q and 77 (in A) variables 
are restricted by the two theta functions appearing in the integrals in Eqs. (|2fj|-p7|). The integrals involving V{m and 
V 1 S 112 are performed within the planes formed by those regions in the q space where 

0[-A]-0[-£i(k + q)] ^0; (30) 

while the integrals involving V^ 1221 and are calculated within the planes defined by 

f[-A- Wo |-SK 2 (k + q)]^0. (31) 

In the same way, the integrals involving V 2222 and V 2221 are performed within the planes defined by 

0[-A]-0[-£ 2 (k + q)]^O; (32) 

and, finally, for V 2 S 112 and V^^the integrating plane is defined by 

9 [—A + lj ]- 6 (k + q)] ? 0. (33) 

The inelastic scattering rates in the Eqs. (pC ■ 27 ) vanish outside each corresponding integrating plane, which means 
that the momentum and energy conservation cannot be simultaneously obeyed for such values of (k, k + q) , and 
therefore no Coulomb scattering is allowed there. It is easy to see that these integrals are non-vanishing only if the 
corresponding integrating plane contains either some part of the single-particle excitation continuum or some branch 
representing the collective excitations (plasmons) in the 2-D q plane. This is of course expected since a finite scattering 
rate must involve real excitations, which in this case are single-particle and collective plasmon excitations. 
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III. NUMERICAL RESULTS AND DISCUSSIONS 



A. Coulomb Coupled Bilayers With Interwell Tunneling 

We consider first two coupled symmetric identical quantum wells of same width W\ = W2 = 150 A with an interwell 
tunneling induced by a barrier of height 228 meV and width 30 A. The total electron density N e = n\ + ?i2 = 2 x 10 11 
cm" 2 in all structures studied in this paper, with ri\ and n-i being the density in the subband 1 and 2, respectively. For 
these sample parameters, the Fermi wavevectors in the first and second subband are kp\ ~ 0.88 xfO 6 cm -1 and k S p 2 
~ 0.69 xfO 6 cm -1 , respectively (the superscript sy stands for symmetric). Here, both subbands f and 2 (symmetric 
and antisymmetric, respectively) are occupied by carriers with n\ ~ f .23 x fO 11 cm -2 and 712 — 0.77 x fO 11 cm" 2 . 

The plasmon dispersion relation is determined by the roots of the dctcrminantal equation det \siji m (q, lo)\ = 0, 
which, after some algebra, can be rewritten as 

£intra£inter ~ [(1 " VaU^V^ + (l - VbII° 2 ) V^U° n - 



2VcVjV H U° n Il 22 (n? 2 + U° 21 )] = 0, (34) 

where 

Smtra = (l - K^) (l - V B H° 22 ) - V^^ (35) 

and 

Sinter = I - V D (ll 12 + H° 21 ) . (36) 

For notational simplicity, we do not explicitly write the energy and wave vector dependence in Eqs. (|34|-j36"l) . For the 
present symmetric situation, the unscreened Coulomb potential Vj = Vh = by virtue of parity symmetry, because 
the wavefunctions (j>i(z) and <fe(z) & r e symmetric and antisymmetric functions of z, respectively. According to the 
Eq. (|34|), therefore, the plasmons dispersion relation in our symmetric bilayer structure are determined by the roots 
of the equation £intra£inter = 0, i.e. either £i n tra = corresponding to the 2D intrasubband plasmons, or Sinter — 
corresponding to the intersubband plasmons. 

There are four roots of Emtra = 0. Two of them are shown in Fig. 1(a) by the solid lines indicating the intrasubband 
plasmon modes (1, 1) and (2,2). Notice that, for each solid line, there is a corresponding dashed line which is also 
the root of the same equation always lying in the corresponding single-particle excitation continuum. It is well known 
that the plasmon modes indicated by the dashed lines inside the single-particle continua are strongly Landau damped 
by single-particle excitations and will be ignored in the following discussion. Furthermore, the intersubband plasmon 
mode (1,2) comes from the roots of £i n ter — 0. The wavefunctions 4>i(z) and 4>2{z) are schematically shown in the 
inset by the solid and dot lines, respectively. Notice that, one is always able to separate the intra- and inter-subband 
plasmon modes in structures which are invariant under space inversion. In addition, the intrasubband plasmons are 
not Landau damped by intersubband single-particle excitations and vice versa in symmetric bilayer systems. The 
single-particle continua SPEn (intrasubband SPE) and SPE12 (intersubband SPE) in Fig. 1(a) are those regions 
where Im {ll] ) 1 (g, cj)} 7^ and Im {Hl 2 (q,w)} 7^ 0, respectively. For the sake of simplicity, we will not indicate 
the continuum SPE22 in this paper because it lies totally inside the continuum SPEn. Moreover, we claim that the 
plasmon mode (2, 2) should be strongly damped by single-particle excitations in the SPEn continuum and will also be 
ignored in the following qualitative scattering rate discussion. Our numerical results of course include all contributions 
as obtained by evaluating the 2D integrals in Eqs. ( f2"of|27j ) . 

Fig. 1(b) shows the same plasmon dispersion relation as in Fig. 1(a) but now in a two coupled asymmetric 
quantum wells with an interwell tunneling. Here, one well is of width 150 A and the other is of width 140 A. For these 
parameters, the Fermi wave vector in the subband 1 and 2 are k F1 ~ 0.93 xlO 6 cm" 1 and k F2 ~ 0.62 xlO 6 cm" 1 , 
respectively (the superscript a stands for asymmetric). Both subbands are occupied with ri\ ~ 1.37 xlO 11 cm -2 and 
ri2 — 0.63 xlO 11 cm" 2 . In this asymmetric situation, the plasmon modes are obtained directly from the roots of 
Eq. d34|). We show in Fig. 1(b) all these roots. We mention that it does not make sense naming the solid lines as 
pure intra- or inter-subband plasmon modes because the structural asymmetry leads to a strong coupling (or mixing) 
between them, and this intra- inter mode coupling eliminates the simplicity of Fig. 1(a). The solid line that is of 
finite frequency as q — > in Fig. 1(b) is the intersubband-like plasmon mode (1,2). This mode enters the continuum 
SPE12 at q ~ 0.42 xlO 6 cm" 1 and should be, in principle, Landau damped. For small values of q, we find the same 
number of roots as in the symmetric situation. The interactions Vj and Vh are finite in the asymmetric case and are 
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responsible for the strong mixing between the intrasubband-like plasmon mode (1,1) and the intersubband-like mode 
(1,2) around q = 0.18xl0 6 cm -1 . Moreover, when the asymmetry is introduced, the depolarization shift (i.e. the 
shift of the intersubband plasmon from the subband ener gy s eparation E 2 i) in the intersubband-like plasmon (1,2) 



at q = decreases. We point out that these roots of Eq. (33) do not provide a complete description of the plasmon 



modes in asymmetric bilayer structures. A detailed theoretical calculation of the dynamical structure factor giving 
the plasmon spectral weight provides a complete picture of the collective mode spectra and can be obtained using our 
multisubband theory. 

Having studied the plasmon dispersion relations we now investigate in Fig. 2(a) the corresponding total inelastic 
Coulomb scattering rate ai(fc) (thick-solid line) and 02(fe) (thick-dashed line) of fast electrons in the subband 1 and 2, 
respectively, as a function of wavevector k in our symmetric bilayer structure. The symbols on the thin lines identify 
the contributions to o\ (k) and (72(h) coming from the emission of single-particle and collective excitations individually. 
The dynamically screened Coulomb interaction components entering in Eqs. j20|-p7|) can be calculated from Eq. ( |T3] ) 
and (|l4|). After some algebra, we get 

_ v A (i- v B n° 22 ) + v§n° 22 

V im — - , {£{) 

£ intra 

_ y- B (i-y- A n° 1 ) + ygn; 1 
^2222 > w°) 

tintra 

and 

V D 

^1221 = ^2112 = ~ • (39) 

Winter 

For the symmetric well case the off-diagonal components of the Coulomb potential all vanish by parity: Vf 112 — 
V n21 = V1211 = V 01 11 = V2221 = V2212 = ^2122 = ^1222 = because Vj = Vr — for symmetric systems. Therefore, 
according to Eqs.( po[p9| ), the total inelastic scattering rates in the subband 1 and 2 are 

ai(k) = ann + C1221 (40) 



and 



<7 2 (fc) = (72222 + CT2112, (41) 



respectively. The terms aim, C122I) 0-2222 and 0-2112 involve integrations of the interactions VA ; , V{ 221 , Villi and 
V 2112 , respectively; in Eqs. (pOj), (|2l|), ( p7| ) and (p6|). The self-energy components in Eqs. d22]-p5|) are zero in the sym- 
metric case. Intrasubband contributions to the scattering rates arise from the terms aim and 02222 , while intersubband 
contributions are due to the terms ai22i and 0-2112 ■ The contributions coming from the plasmon modes (filled-square 
lines) are obtained separately by excluding the continua SPEn and SPE12 from the numerical integrations, whereas 
contributions coming only from the single-particle continua are obtained by numerically evaluating Eqs. ( ^l| ) and ( pfij ) 
only for the region representing each continuum leaving out the plasmon contributions. Single-particle excitations 
contribute for all values of wave vectors k. However, neither intra- nor inter-subband plasmon mode contributes to 
the scattering rates close to fc^j or k S p 2 . These collective modes provide excitation channels for inelastic relaxation 
only above some threshold wavevectors. The intrasubband plasmon mode (1, 1) begins to contribute to either ai(fc) 
or a 2 (k) when the wavevector is larger than the same threshold ~ 1.65xl0 6 cm -1 . On the other hand, the 
contribution coming from the plasmon mode (1,2) has a different threshold for each scattering rate. This mode 
begins contributing to ai(fc) and a 2 (k) when the wavevector is larger than the thresholds k^ 2 ~ 1.25xl0 6 cm -1 and 
kl 2 , ~ 2.0 xlO 6 cm -1 , respectively (notice that k S p 2 < k S p Y < fc*| < ^11 < ^12')- Obviously, these thresholds depend 
on the particular choice of sample parameters. In the present paper, they are smoothed ( instead of being a very 
sharp threshold ) because we are considering the impurity induced constant 7 = 0.2 meV in our numerical evaluation. 
These thresholds become much sharper for smaller values of 7 without any other substantive changes in our numerical 
results. 

Figure 2(b) shows the same results as in Fig. 2(a) but for the asymmetric bilayer system of Fig. 1(b). In contrast to 
the symmetric case, where we were able to separately obtain the inter- and intra-subband plasmon modes through the 
roots of Sinter = and Sintra = 0, respectively; the coupled plasmon dispersion in the asymmetric system is obtained 



directly from the numerical roots of Eq. ( 34 ) in which the bare off-diagonal Coulomb interactions Vj and Vh are now 
non- vanishing. The terms in Eq. ( |34| ) involving Vj and Vh are responsible for the mixing between the inter and 
intra-subband plasmon modes and for not allowing the contributions coming from the intra- and intersubband-like 
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plasmon modes (1, 1) and (1,2) to be picked up completely separated from each other in the scattering rate. Notice 
that dynamically screened Coulomb potential V^ lm is a full 16x16 matrix (for the 2-subband model — in general, 
it is an n 4 x n 4 matrix for an n-subband problem) in the present situation and is obtained from Eq. (|l3|), which 
involves the dielectric matrix Siji m {q, to) and the bare Coulomb interactions Va,Vb, Vc, Vd, Vj, and Vh (all of which 
are finite in this strongly coupled asymmetric bilayer structure) M. Therefore, both inelastic scattering rates ai(k) 
and <72 (k) in the asymmetric case contain all terms shown in Eqs. ( f20f]27| ), which are finite in this situation. For the 
sake of clarity and to understand Fig. 2(b) in the same way we did the details in Fig. 2(a), we choose to show three 
contributions to the inelastic scattering rates a\{k) and o~2{k) in Fig. 2(b) separately: the single-particle excitations 
in the continua (i) SPE12 (triangles-up) and (ii) SPEn (diamonds); and (iii) the plasmon mode segment outside these 
continua (filled squares). The filled squares in Fig. 2(b) represent contributions coming from those segments of the 
plasmon modes which lie outside any single-particle excitation continua (see Fig. 1(b)). Contributions coming from 
the plasmon segments lying inside each continuum have been kept in our numerical work along with the single-particle 
excitation contributions because it is essentially numerically impossible to separate the two in this regime. We should 
mention that, due to the fact that one is not able to eliminate the contributions coming from the over-damped plasmon 
modes lying inside the Landau continua, the thin lines in Fig. 2(b) serve only the purpose of qualitative illustration 




B. Coulomb Coupled Bilayers With No Interwell Tunneling 

Now we investigate the two Coulomb coupled identical symmetric quantum wells of width W\ = W2 =150 A each 
with no interwell tunneling (i.e. the interwell barrier is taken to be infinity) and with a barrier width of 30 A. Here, 
the Fermi wavevectors in the two wells are of the same value, i.e. fepo = &fi = kp2 — 0.79 x 10 6 cm -1 (or, equivalently 
n>i — TI2 — 10 11 cm -2 ). Notice that, the indices 1 and 2 should now be treated as well indices since there is no 
tunneling induced bonding-antibonding states. As we mentioned in the Introduction, an energy degeneracy arises 
in this case, i.e. E\ = E 2 because the two wells are identical with no interwell tunneling. If there is no tunneling, 
the bare Coulomb potential components Vj = Vh = Vd = and the polarizability II^ = = independent of 
whether the bilayer structure is symmetric or not. Besides, for this symmetric no-tunneling bilayer structure, the bare 
Coulomb potential Va = Vb by symmetry and the polarizability VP 1X = II^ = IIo due to the fact that the densities 
in each well are identical. 

According to the Eqs. ( p4f|36| ) , therefore, the plasmon dispersion relation should be obtained only from the roots of 

e% = (1 - V A I1 ) 2 - V^Ul = 0. (42) 

Here, the subscript (superscript) nt (sy) stands for no tunneling (symmetric). As shown in Fig. 3(a), we find four 
roots of Eq. (fL?|). The solid curves correspond to the in-phase optical, uj + (q), and the out-of-phase acoustic, ui-(q), 
plasmon modes in the bilayer structure pd| ] . These u>± (q) modes have been observed p"2[ in multilayer semiconductor 
systems via inelastic light scattering spectroscopic experiments. They represent in-phase and out-of-phase interlayer 
density fluctuation modes: the out-of-phase acoustic mode, to~{q — > 0) ~ 0{q) represent densities in the two layers 
fluctuating out of phase with a linear wave vector dispersion and the in-phase optical mode, uj + (q — ► 0) ~ y/N e q, 
represent densities in the two layers fluctuating in phase with the usual 2D plasma dispersion. The dashed lines 
represent those collective modes which should be strongly Landau damped by the single-particle excitation continuum 
SPE, i.e. the region where Im {Ho(g, uj)} 7^ 0. 

Fig. 3(b) shows the same as in Fig. 3(a) but for an asymmetric no-tunneling situation with the two wells being 
different, one with a width of W\ =150 A and the other a width of W2 =142.4 A. In this case, our no-tunneling 
bilayer structure is no longer invariant under space inversion and, consequently, the energy level degeneracy is broken, 
leading to the energy E\ < 2?2- Besides, the bare Coulomb potential Va is now different from Vb- As we discussed 
before, the two wells now have different charge densities but we consider the whole system still being in equilibrium. 
Furthermore, the Fermi wavevector in the first and second subband are the same as indicated before, i.e. k F1 and k F2 , 
respectively. Because of the densities in the two wells being different from each other, the polarizability H^ ^ n^. 
The shadow area in Fig. 3(b) is the single-particle excitation continuum in the wider quantum well, i.e. the region 
where Im {n^g, cj)} 7^ 0. The plasma dispersion relation is now given by the roots of the Eq. (p5|). Note that 
all plasma modes in the zero tunneling system are by definition intrasubband plasmons in our model where higher 
subbands are neglected. 

As shown in Fig. 3(b), we again find four roots of such an equation and consider that the dashed lines should be 
strongly Landau damped modes since they are inside the single particle continua. Furthermore, it does not make sense, 
in principle, to define the solid lines in Fig. 3(b) as pure acoustic or optical plasmon modes because the asymmetry 
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leads to a difference between the electron densities in each layer. Now, the wider well has 30% more electrons than 
the narrower one and, consequently, the densities in the two layers are not fluctuating exactly either in phase or out 
of phase. The solid lines in Fig. 3(b) are the approximate acoustic- and optical- like plasmon modes with the strict 
distinction meaningful only in the long wavelength limit. Due to the structural asymmetry the acoustic-like plasmon 
mode enters the SPE continuum at a smaller wavevector leading to significant Landau damping of the acoustic plasmon 
mode by single-particle excitations in the asymmetric bilayer system. In the single-particle continuum of the layer 2 
(the narrower well) the acoustic-like plasmon mode is completely suppressed and we find no acoustic-like mode in the 
Im{ri2 2 (<7, to)} 7^ regime. In general the acoustic-like plasmon mode is found to be much more sensitive to small 
asymmetry effects than the optical- like plasmon mode p3| . This is physically reasonable and should be experimentally 
tested via inelastic light scattering experiments. 

Now we concentrate on the investigation of the scattering rates eri(fc) and cr 2 (fc) in the symmetric bilayer structure 
with no tunneling. As the bare Coulomb potential Vj ~ Vu = Vd = 0, it is straightforward to see that only Vf xll 
and V2222 are finite in the screened Coulomb interaction matrix V s for a bilayer structure without any tunneling. 
Therefore the scattering rates in Eqs. ( plf|2g) all vanish by symmetry in this case. The only nonvanishing terms to 
be calculated are aim and 02222 in Eqs. ( ppj ) and (p7|), respectively. Furthermore, as we discussed before, we have 
the polarizability 11^ = LI" , = lip and the bare potential Va = Vb for identical (i.e. symmetric case) quantum wells. 



According to the Eqs. ( |37|) and (38), therefore, the screened Coulomb potential is given by 



V A (1 - v A n°) + v2n° , x 

V?1H - Kf 222 - ^sy 2 — (43) 

in the present situation of a symmetric bilayer system with no intcrwcll tunneling. In fact, the total inelastic Coulomb 
scattering cri(k) and (J2{k) are identical because the two wells are identical with the same density. The thick line 
shown in Fig. 4(a) represents the total inelastic scattering rate which is equal (ci(fc) = cr ml = ^(fc) = 02222 ) m 
both subbands, as a function of the wavevector k. To show separately the contributions coming from the emission of 
plasmons (squares) and single-particle excitations (diamonds), we again exclude the region where Im [Ilo((?, u>)] =/= 
from the numerical calculations to obtain the plasmon contribution. Single-particle excitations again contribute at all 
values of the wavevector whereas the plasmons begin contributing to the scattering rate for wavevectors k larger than 
the Fermi wavevector kpQ. There are clearly two thresholds wave vectors in the plasmon contribution (squares), one 
at k = ~ 1.25x 10 6 cm -1 and other at k = ~ 1.65x 10 6 cm -1 . These are the thresholds for the emission of the 
acoustic and the optical plasmon, respectively. The substantial difference between Figs. 4(a) and 2(a) demonstrates 
the strong effect of tunneling on the inelastic scattering rates in bilayer structures. This is one of the important new 
qualitative results in our paper. 

Figure 4(b) shows the same results as in Fig. 4(a) but for the asymmetric bilayer structure without tunneling. 
As we discussed before, the asymmetry leads to II^ 7^ n^. Furthermore, the bare Coulomb potential Va 7^ Vg 
and, therefore according to Eqs. (|3^) and fl3S|), the screened Coulomb potential Vf ln 7^ V2222 m tne asymmetric 
case. In this situation, ci(fc) = trim (thick solid line) and U2{k) — 02121 (thick dashed line) represent the total 
inelastic Coulomb scattering rates in the wider and narrower layer, respectively. They are different from each other 
because the two wells have different widths and densities in the asymmetric situation. Again, we separate the different 
contributions (by plasmons and by SPE) by excluding the single-particle excitation continuum SPE from the numerical 
calculations to obtain the plasmon contribution. It is important to point out again that the squares in Fig. 4(b) 
represent contributions coming from the emission of undamped plasmon modes whose frequency uj(q) lies outside the 
continuum SPE (see Fig. 3(b)). There is only one threshold wavevector k ~ 1.71 xlO 6 cm -1 in the thin solid line 
(squares) corresponding the plasmon contribution to oi(fc). This threshold is due to the emission of the optical-like 
plasmon mode shown in Fig. 3(b). We also find that, the thin solid line (diamonds) corresponding to the SPE 
contribution to ai(k) does not contain any contribution coming from the acoustic-like plasmon mode. As a matter 
of fact, there is no contribution to a\{k) in Fig. 4(b) coming from the emission of the acoustic-like plasmon mode at 
all because the integral in u\{k) does not contain any segment representing the acoustic-like plasmon mode which is 
heavily Landau damped in the asymmetric situation under consideration. On the other hand, the thin dashed line 
(squares), corresponding to the plasmon contribution to &i(k), clearly has two threshold wave vectors k ~ 1.15xl0 6 
cm -1 and k ~ 1.76xl0 6 cm -1 , which characterizes the emission of the acoustic- and optical-like plasmon mode, 
respectively. Thus, in the asymmetric case, the acoustic-like plasmon modes contribute to carrier scattering 02 (A;) in 
the narrower well but not to <7\{k) in the wider well by virtue of strong Landau damping. The difference between 
Figs. 2(b) and 4(b) represents the strong effect of tunneling on the second component of the inelastic scattering rates 
cr 2 (&;) in bilayer asymmetric structures. 



C. Single Symmetric Quantum Well 
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We now consider (for the sake of comparison) a single symmetric GaAs-Al x Gai_ x As quantum well of width 300 
A, barrier height 228 meV, and with the same total electron density N e = 2xl0 n cm~ 2 as used before. These 
sample parameters lead to the Fermi wave vector in the first subband A;™ 8 ' ~ 1.13x 10 6 cm _1 with only one subband 
occupancy. Here, the second subband is empty, which leads to I±2 2 = 0. As we discussed before, only the bare 
Coulomb potential Va,Vb , Vc and Vd are finite because Vj = Vh = in symmetric structure. According to Eqs. 
([35|) and (39), therefore, the intra- and inter-subband plasmon modes are obtained from the roots of the equations 

4t: = (1-^0 = 

and 

43£ = i-Mn?2 + n8i)=o, 

Taking II^ = (unoccupied excited subband) in Eqs. (p57l), (p^) and (59h we get 



^111 = ( 44 ) 



"intra 



_ Mi-^n? ) + ^n; 

^single v ' 



intra 



and 



^221 = ^112 = =• (46) 



Again in this case, the screened Coulomb potential Vf 112 = Vf 121 — Vf 211 — V 2 \u = V 2221 = ^2212 = ^2122 = ^1222 = 
by symmetry because Vj = Vh = 0. Therefore, as we discussed in the Sec. III. A, the total inelastic scattering rates 



in the first and the second subband are given by Eqs. (|40j) and (41), respectively. In the same way we did for the 
bilayer structures, we present the scattering rates 01 (fc) and 02 (A; ) in Fig. (5) in the thick-solid and thick-dashed lines, 
respectively. The symmetric nature of the single well system enables us to separate out the different contributions to 
the scattering rates as discussed before. We find that contributions to eri(fc) come mainly from the emission of both the 
intrasubband plasmons (1,1) and the intrasubband single-particle excitations SPEn. The emission of intersubband 
excitations turn out to make negligible contributions to the scattering because of the sufficiently large energy gap 
between the two subbands (loq = E 2 \ = 14.63 meV). For this particular choice of the sample parameters, axm turns 
out to be much larger than (71221 , implying that the carrier relaxation process in the ground subband is almost entirely 
via intrasubband scattering. Another important point in Fig. 5 is that inter- and intra-subband plasmon modes as 
well as intra- and inter-subband single-particle excitations contribute to the total inelastic scattering rate a 2 {k) in 
the second subband. Notice that, in contrast to the behavior of ci(fc), the total inelastic scattering rate in the second 
subband <r 2 (k) does not vanish for wavevectors any k. This is due the fact that there is no Fermi surface in the second 
subband . This should lead to qualitatively different effects in the measured carrier injected in the second subband 
compared with that in the ground subband fl4|| . This lifetime, which is inversely proportional to the total inelastic 
scattering rate a 2 {k) 1 should be finite for all finite wavevectors in the excited empty subband. 



IV. CONCLUSIONS 



We have developed a theory for calculating the inelastic relaxation rate for Coulomb scattering in coupled bilayer 
structures in semiconductor double quantum well systems. We use a many body theory based on a multisubband 
generalized GW approximation with the inelastic scattering rate defined by the magnitude of the imaginary part of 
the on-shell electron self-energy. Effects of dynamical screening, mode coupling, and Fermi statistics are naturally 
included in our many-body theory. We demonstrate the usefulness of our theory by obtaining results for general 
representative two-subband model systems: Coulomb coupled bilayer GaAs-AlGaAs double quantum well structures 
both with and without interwell tunneling and also with and without interwell asymmetry in the system. Our theory 
naturally allows for distinguishing various physical mechanisms contributing to the inelastic scattering rate: Intra- and 
inter-subband contributions. We provide a critical qualitative discussion of these various contributions to scattering 
and comment on the effect of interwell tunneling and structural asymmetry in bilayer quantum wells. 
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FIG. 1. Plasmon dispersions in two coupled GaAs/Alo.3Gao.7As quantum wells of widhts (a) Wi = W2 = 150 A (symmetric) 
and (b) Wi = 150 A and W2 = 140 A (asymmetric); and separated by a barrier of width 30 A . For the symmetric (asymmetric) 
situation the energy separation between the two subbands is luq — 1.68 meV (u>o = 2.62 meV). The shadow areas present the 
single-particle excitation regions SPE n „/ where Im{n° n ; (q, w)} 7^ 0. Each structure is shown in the inset where <f>i(z) and 
c(>2 (z) are schematically shown by the solid and dot lines, respectively. 

FIG. 2. Total inelastic Coulomb scattering rate of electrons in our coupled bilayer (a) symmetric and (b) asymmetric 
structures. The thick-solid and thick-dashed lines denote the total scattering rate a n (k) for n=l and 2, respectively. The 
symbols represent the different contributions to the scattering. 

FIG. 3. Plasmon dispersions in two coupled with no interwell tunneling GaAs/Alo.3Gao.7As quantum wells of widths (a) 
Wi — W2 = 150 A and (b) Wi = 150 A (symmetric) and W2 = 142.4 (asymmetric) A ; and separated by an infinity barrier 
of width 28.87 A . The shadow areas present the single-particle excitation regions SPE n / where Imln^/ (q, uj)} 7^ 0. Each 
structure is shown in the inset where (fri(z) and <j>2{z) are schematically shown by the solid and dot lines, respectively. 



FIG. 4. Total inelastic Coulomb scattering rate of electrons in our coupled bilayer (a) symmetric and (b) asymmetric 
structures with no tunneling. In the part (b), the thick-solid and thick-dashed lines denote the total scattering rate a n (k) for 
n=l and 2, respectively. The symbols represent the different contributions to the scattering. 
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FIG. 5. Total inelastic Coulomb scattering rate of electrons in a single quantum well. The thick-solid and thick-dashed 
lines denote the total scattering rate a n (k) for n=l and 2, respectively. The symbols represent the different contributions to 
the scattering. 
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